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^ ■ Abstract 

We have calculated the effect of the inclusion of final state interactions on pion 
production from nuclei with incident pion beams. We find that the effect of absorp- 
tion, along with geometric considerations, explains much of the enhancement at low 
I invariant mass seen in recent data for invariant mass spectra with increasing atomic 

QQ ■ number. While the variation among nuclei is well reproduced for all charge states, the 

. enhancement for the ratio for heavy nuclei to deuterium for the 7r"*"7r~ final state is 

^ I only partially understood. 



1 Introduction 



The pion induced pion production on both a free nucleon and a nucleon within nuclear 
material provides an opportunity to study the mechanism for pion production. In a 
series of papers ||l], ^, ^, ^ the results of experiments done at the CHAOS (Canadian 
^ ' High Acceptance Orbit Spectrometer) spectrometer at TRIUMF were presented on 

■ pion production with pion beams for several nuclei. Specifically, the collaboration 

reported data for the n^A — > tt^tt^X and tt^A — > vr+vr+A reactions on ^H, ^^C, ^^Ca 
and ^''^Pb at an incident pion kinetic energy of 283 MeV. The results showed that the 
shape of the invariant mass distribution of the two final pious changes with atomic 
mass number. 

An examination of the changes among the heavy nuclei (carbon and heavier) reveals 
that the variation is about the same for the two final charge states. The shape of the 
spectrum for the case of a deuteron target is very different for the two cases. 

A striking feature of the raw spectra is a large peak at small invariant mass of the 
two-pion system. This peak is at least partly the result of the experimental conditions 
and can be regarded as a Jacobian peak related to the acceptance of the spectrometer. 
Since the correction for this large instrumental effect is difficult to make, the results 
are often presented as ratios, in the expectation that the acceptance effect will cancel. 
Thus, Ref. @ gives ratios to the deuteron cross section, a common way to normalize. 
However, in this case, since the deuteron results are very different for the two final 
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charge states, this procedure results in very different spectral ratios for them. If one 
instead takes ratios to the carbon cross sections, the calcium and lead results are about 
the same for the two final charge states. 

Shortly after these first data, experiments Q done with the Crystal Ball at Brookhaven 
National Laboratory produced data on the production of two neutral pions for D, C, 
Al, and Cu with a negative pion beam. One sees a very similar trend for the ratios in 
these data to that of the CHAOS data. 

Since the basic conclusion of these data sets is that the heavier the nucleus, the 
more the invariant mass spectrum is peaked toward lower values of invariant mass a 
number of groups have seen in this data possible evidence for a supposed modification 
of the TTTT interaction in nuclear matter. 

In early work, a theoretical expectation of an enhancement of strength in the vr+vr^ 
channel for I=J=0 near the threshold was suggested by Schuck et a/.0] based on an 
analogy with Cooper pairs. 

Recently, several studies have placed emphasis on medium effects to explain the 
strength enhancement of M^t^ distribution in vr+vr^ channel {a -meson channel). The 
calculation of Rapp et al. [Q] is based on a simple model in which the pion production 
takes place as an elementary reaction. In their model two contributions were included, a 
single pion exchange reaction and the excitation followed by the decay of the A^*(1440) 
into two pions and a nucleon. The vrvr final state interaction is included by using the 
chirally improved Jiilich modelQ and several medium effects were included. There is 
reasonable agreement of their calculation and the CHAOS data in both the vr+vr" and 
Tr~^TT~^ channels. 

Vicente- Vacas and Oset used a macroscopic model for the pion production with 
final state interaction among the pions included in the nuclear medium. Several nuclear 
effects were taken into account in their study including an approximate treatment 
of pion absorption, Pauli blocking and Fermi motion. Their calculation reproduced 
the CHAOS data for calcium in the vr+vr"'' channel reasonably well and also gave a 
reasonable agreement for the deuteron. In contrast, this model failed to reproduce the 
CHAOS data in 7r"'"7r~ channel except for ^H. 

The enhancement of strength in the I=J=0 channel vr+vr" channel near the thresh- 
old was argued by Hatsuda et al.|ll] to be due to the partial restoration of chiral 
symmetry in nuclear matter. 

Davesne et g/.p^j included the effect of chiral symmetry restoration and the in- 
fluence of collective nuclear pionic modes and found an enhancement of the spectral 
function a small energy. 

Aouissat et a/.|13] argued that the combination of effects of the restoration of chiral 
symmetry |11| in nuclear matter and standard many body correlation could explain 
the strength of the spectrum M^t^ near threshold |^. 

While these considerations may be correct there may also be more prosaic reasons 
for the increasing strength at low invariant mass. We investigate the effect of the flnal 
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state interaction of the two individual pions with the nucleus, especially the absorption 
of the pion. Our calculation contains no free parameters and pays special attention to 
the geometry of the reaction. 

We make two simplifying assumptions. 

1) We take the reaction to be coherent on the nucleus. That is, the nucleus is 
assumed to carry off the remaining energy and momentum as a single entity. Since 
it is much heavier than the pions, it can carry a large amount of momentum with 
negligible amount of energy. In this extreme case the total beam energy is passed on 
to the two-pion system. The invariant mass of the two-pion system is given by 



where is the pion mass (charged or neutral as appropriate) and q is the momentum 
of one of the pions in the two-pion rest frame, the internal momentum. The above 
assumption then leads to the relation 



u; = ^Af2 + Q2 (2) 

where uj is the initial beam energy and Q is the total momentum of the two-pion pair. 

2) We assume that the production takes place at a fixed point on the surface of the 
nucleus in a plane at the equator perpendicular to the beam direction. The reaction 
probability is greatly enhanced in this region for two reasons. First, there is a simple 
geometrical factor of the distance to the center of the nucleus. Second, if the incident 
pion comes in with a small impact parameter either the initial pion will be scattered 
(even a moderate loss of energy will make the pion production impossible) or one of 
the final pions will very likely be absorbed, having to pass through the entire nucleus 
to escape. Thus, the reactions in which the incident pion has a small impact parameter 



are greatly suppressed. This point is discussed in Ref. [10| 



2 Technique 

2.1 Qualitative Overview 

It is possible to see in a simple qualitative manner that pion absorption in the final 
state will lead to a relative enhancement of the invariant mass spectrum at low pion 
invariant masses. 

For a first orientation we will assume that the pion pair proceeds forward in the 
direction of the initial beam since the distribution of the momentum of the pair is 
expected to be very forward peaked. For the purposes of the diagram only (see Fig. 
|l]) we also take the momentum of one pion in the frame of the two pion system, q, to 
be perpendicular to the beam direction. 
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With these simphfying assumptions we may now consider the effect of the absorp- 
tion of the pions in the final state on the shape of the invariant mass spectrum. Figure 
|l] shows two cases. The left portion of the sketch corresponds to the situation where 
the invariant mass is small, so the internal momentum is small and the momentum of 
the total two-pion system is large, leading to a narrow cone representing the paths of 
the two pions. It is seen that the pions are exposed to relatively little interaction with 
nuclear matter. In the limit of the internal momentum becoming zero, the pions would 
encounter no higher nuclear density than that in which they were formed. In this limit 
the size of the nucleus would matter very little. 

In the right portion of the diagram a case is illustrated in which the invariant mass 
and internal momentum are larger. In this case it is seen that a large fraction of the 
pions must traverse a significant amount of nuclear matter. The larger the nucleus, 
the larger the probability that one member of the pion pair would get absorbed. 

Hence, we can anticipate that the geometry of the absorption in the final state leads 
to a depletion of the higher values of the invariant-mass spectrum and that the effect 
becomes more important with increasing A. 



2.2 Calculation 

The method used here is very similar to that of Oset and Vicente Vacas [14]. However, 



they calculated the effect only for calcium and used an eikonal approximation for the 
propagation of the final pions. 

We assume that the pion pair is formed quickly at short range by the strong in- 
teraction, and that this formation is independent of the nucleus that it is near. Also, 
we assume that the final total momentum of the two pions is in the same direction 
as the incident beam. This assumption of exact forward propagation of the pair was 
investigated by calculating with finite angles. It was found that the dependence on 
this angle was very small. We define the spherical angles 6 and (p from the beam axis 
as the direction of the tt"^ in the center of the mass of the pair. 

After formation the pions are allowed to propagate independently in Coulomb and 
strong nuclear potentials with relativistic motion. 

Their momenta and coordinates are propagated with the equations 

p(t + 6t) = p{t) + FSt (3) 

and 

r(t + 6t) = r{t) + lll^i (4) 

where F is the resultant of the Coulomb and strong interaction forces. We take the 
electric charge distribution of the nucleus to be 

Pc{r) = (5) 

1 + e a 
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The values of 7 and a were taken from Ref. |jT^ and are (in fm) (2.44, 0.5), (3.07, 
0.519), (3.65, 0.54), (4.252, 0.589) and (6.6, 0.55) for C, Al, Ca, Cu and Pb respectively. 
The constant po was fixed by normalizing the integral of the density to the total charge 
Ze. 

The nuclear potential is given by 

27r 

Vs{r) = T (iV/.„ + Zf^p) p{r) (6) 

mA 

where N, Z are the number of neutrons and protons and A is the atomic mass number. 
p{r) is the nuclear strong density, i.e. the density of the centers of the nucleons. 

p{r) = (7) 
1 + e~ 

The values of c and a were (2.25, 0.5), (2.95, 0.5), (3.55, 0.54), (4.23, 0.55) and (6.5, 
0.55) for C, Al, Ca, Cu and Pb respectively and f-j^n and /vrp are the real parts of the 
forward amplitudes of the vrn and vrp amplitudes [p^. This real potential is modeled 
after the pion-nucleus optical potential. We will see later that these non-absorptive 
strong interactions play a minor role in the result. 

For each initial value of the invariant mass {Mi) we calculate the internal momen- 
tum, the total momentum of the two-pion pair and, with values of the angles of these 
two quantities, the initial values of the momenta for each pion can be determined. The 
pions are then propagated to a large distance (into free space) and a new invariant 
mass is calculated for the pair (Mf). In this way the function 

Mf{x,<P,Mi), (8) 

where x = cos 6, is established. It depends on the angles of both q and Q but we 
suppress the dependence on the angles of Q. 

The final spectrum (without absorption) will be given by 

Because the derivative in the denominator is sometimes zero, sharp peaks appear in 
the spectra for individual values of the momentum directions. The integration over 
the angles results in a smooth spectrum but, since the integration is done numerically, 
some care is required to obtain a presentable curve. The lines shown in the figures 
have been smoothed since the final calculated values still show a small residual ripple. 
We use as a simple model for the initial spectrum two body phase space. 



2.3 Absorption Factor 

Pion absorption is an important factor during the passage of pions through the residual 
nucleus. It reduces the cross section for pion production and affects the shape of 
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spectrum. We include the effects of absorption by calculating the factor the probability 
that the pion has not been absorbed, P(r), at the point r as follows. The probability 
of survival is given by : 



P(r) = 1 - / X{r')P{r')dl (10) 



where fq is the formation point, and the initial value of P = 1. The integral is taken 
along the trajectory of the particle. As in Ref. the absorption probability per 

unit length is assumed to be proportional to the square of the nuclear density 

A(r) = XoiE)p\r), (11) 

The Ao parameter depends on the energy of the pion and is obtained by fitting to 
absorption data. For light nuclei we first calculated the eikonal expression for the total 
absorption cross section 

a(Ao) = 27r ada {l - e" ^-^ ^^^op^Cv^^) j _ (-^2) 



By comparing the values of o"(Ao) with experimental values of (Jabs{E) [O, 18, 19 



21, 22 1 we inferred the values of Xo{E) needed to fit the data (see figure |2|). We found 



that the values of Xo{E) had little dependence on the atomic mass number A for light 
nuclei. 

For heavy nuclei it is important to follow the trajectory instead of using Eq. 
Figure ^ shows the result of the fit for gold. The Aq used in this case is also given in 
Fig. I 

Denoting the result of following the trajectory of the i*^ pion into free space by 
Pi{x,(p,Mi) the final spectrum is given by the integral over all possible decay angles, 

I rl I-2TV 
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Sf{Mf) = — I dx I d^Sf{x,(^,Mf)Pi{x,<P,M,)P2{x,^,M,). (13) 



3 Results 

Figure ^ shows the resulting spectra as well as the original model spectrum (solid 
line). Also shown is the spectrum with only the strong and Coulomb potentials. These 
spectra cannot easily be compared directly to the data for two reasons. 

1) Our initial model spectrum represents only roughly the free spectrum for vr+Tr^ and 
not well at all for it~^tt~ . In fact, the free spectrum is not well known. Even if it were, 
we cannot be sure that it is the appropriate starting spectrum in the nuclear case. 

2) There is a large instrumental acceptance correction to be made. 

For these reasons we, like the experimental group, compare with ratios. Since we 
are interested in the increase of the spectrum ratios near small values of invariant mass 
we normalize the ratios (theoretical and experimental) to unity at 330 MeV. 
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The CHAOS group calculated the ratio for both tt+tt and vr+vr"'" to show the 
effect of the nucleus. The strong enhancement near threshold is due to the fact that 
the invariant mass distribution of D (^H) and A (^^C, ^''Ca,^'^^Pb) change dramatically 
near the threshold where is a much smaller peak in the D distribution. 

In our study we calculated the ratio of heavier nuclei to carbon instead of deuterium 
and found that the ratio increases with A for all charge channels by almost the same 
amount. 

Figures |5| and ^ show the ratios of the cross sections for ^''^Pb and ^'^Ca to that of 
^^C. The calculated ratios are seen to reproduce the experimental data in the vr+vr"'" 
and 7r"'"7r~ channels. The two charge channels show similar behavior. Figure |6| also 
shows the effect of varying the assumed radius at which the production takes place. 

The upper part of Fig. ^ shows the difference seen in the lead to carbon ratio 
using the value of Aq determined from the gold data and that obtained from the light 
element fit. The lower part of the same figure shows the difference between using the 
Woods-Saxon density and a shell model density for calcium |23]. 

Figure ^ shows the dependence on the final charge state of the calculations. It is 
seen that the effect is moderate. 

Figure |9| shows the ratio of the cross sections for Cu and Al to C compared to the 
experimental data in the channel. Although the error bars are large, the trend 
of the data is clearly reproduced. 

The estimate of the interaction potential to use in a classical calculation by a 
quantum model must, of course, be approximate. Figure ^ shows that this potential 
plays a minor role by comparing with calculations with it set to zero. 

Figure 11 shows the ratios of the heavy elements to deuterium for the tt'^tt'^ final 
state. The comparison between theory and experiment is satisfactory. 

Figure 12 shows the ratio of the lead spectrum to that of deuterium for the vr+vr^ 
final state. The effect considered here is able to explain only a little less than half of 
the observed experimental effect. 



4 Discussion 

We have seen that the variation in the spectral shape from nucleus to nucleus (above 
deuterium) can be understood in terms of the interaction of the produced pions with 
the nucleus. The size of the nucleus plays a strong role in this explanation. Thus, 
there is little or no evidence in these ratios for a modification of the pion production 
or pion-pion interaction. It is not clear that one should expect an effect of this type 
in any case. As has been discussed, the point at which the production takes place is 
largely determined by the density. Thus, as one moves to heavier nuclei there will be 
a strong tendency for the reaction to occur at the same density. Since the average 
density is the major parameter used to determine the medium effect, and it tends to 
remain constant, one might believe that no nuclear matter effect on the production 
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would be seen in comparing among heavy nuclei. 

Having said that, one must admit that there is some variation in the in the ex- 
pected density going from (say) carbon to lead. If we use the linear dependence of the 
condensate on the density of Cohen, Furnstahl and Griegel [^] with the calculation 
of the change in the shape of the spectrum by Hatsuda and collaborators [0] along 
with a guess of the change in nuclear density, the result is of the same order as the 



discrepancy seen in the upper part of Figs, y or 10. However, one must be very careful 
about trying to draw any conclusion from this since this also about the same order as 
the difference between the two choices in Aq as seen in Fig. |^ and not far outside the 
experimental errors. 

Of course, one might try to compare the production on the nucleon with production 
in the nucleus which is close to what has been done before. One would have to assume 
that the mechanism was simply production from one of the nucleons in the nucleus 
(a reasonable first assumption). In the present experiments the comparison has been 
made with the deuteron. 

It has long been known (see Ref. |22] for example) that the shape of the invariant 
mass spectrum for the vr+vr^ final pair is different from that of the 7r+7r+ pair at low 
incident energy. According to Ref. |1C] the sharp decrease in the low mass region 
is due to a cancellation in diagrammatic contributions. These authors suggest that 
an alteration of the cancellation due to interaction with the surrounding nucleons 
might lead to a modification in this cancellation and thus to a change in shape of the 
production spectrum in the direction seen. 

It is also clear (see again Ref. |2^) that the vt^tt" spectrum shape changes to more 
resemble the vr+vr"'" distribution (and phase space) as the energy increases. As one 
goes from the reaction on a single nucleon to that on a nucleus the effective energy 
at which the reaction takes place becomes higher due to the possibilities of coherent 
production. Thus, one should also expect an increase of strength in the low mass part 
of the spectrum from this effect as well. 

This work was supported by the National Science Foundation (Grant 0099729) and 
the Jordan University of Science and Technology. 
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Figure 1: Schematic view of the basic effect. 
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Figure 2: Aq determined for light nuclei (solid curve) and for gold (dashed curve). 
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Figure 3: Comparison of the absorption cross section computed with the dashed curve in 
figure I with the absorption data of Nakai et al. . 
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Figure 4: Variation of the basic spectrum for carbon, calcium and lead. The solid curve shows 
the shape the initial free spectrum and the long-dashed curve represents the modification 
due to the non-absorptive interactions only. 
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Figure 5: Ratios for lead and calcium to carbon for the tt^tt'^ final state compared with the 
CHAOS data i. 
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Figure 6: Ratios for the tt+tt" final state compared with the CHAOS data (squares) and 
|jl|] (triangles). The solid line is the result of a calculation with r = c and the dashed line 
was made with r = c + a. 
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Figure 7: The upper panel shows a comparison the ratio calculated using values of Aq 
obtained from the fit to gold (solid) and the light nuclei (dashed) . The lower panel compares 
the result using the Woods-Saxon densities given in the text (solid line) with the result of 
using a shell model density for Ca (dashed line). The data is the same as in figure |[ 
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Figure 8: Comparison of results for different cliarge states for the lead to carbon ratio. The 
solid curve corresponds to the tt^tt^ final state the short dash to the neutral final state and 
the long dash to the 7r"'"7r~ final state. 
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Figure 9: Results for the ratios for the n tt final state compared with the Crystal Ball data 
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Figure 10: Results showing the effect of setting the strong potential to zero for the tt"* 
final state. The solid curve is calculated with the strong potential and the dashed curve 
it set to zero. 
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Figure 11: Ratios of the spectra of heavy nuclei divided by the basic spectrum for the 7i~^7i~^ 
final state. The squares and sohd curve correspond to lead, the triangles and long-dashed 
curve to calcium and the diamonds and short-dashed curve to carbon. 
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Figure 12: Ratio of the lead spectrum to the basic spectrum for the 7r~^n final state. 
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